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Abstract 

^ ' The adiabatic approximation is well-known method for effective study of few- 

^I^ . body systems in solid, molecular, atomic and nuclear physics, using the idea of 

^ ■ separation of "fast" and "slow" variables. The generalization of the standard 

dJ . adiabatic ansatz for the case of multi-channel wave function when all variables 

treated dynamically is presented. For this reason we are introducing the step- 
by-step averaging methods in order to eliminate consequently from faster to 

2 \ slower variables. We present a symbolic- numerical algorithm for reduction of 

multistep adiabatic equations, corresponding to the MultiStep Generalization 
> _ of Kantorovich Method, for solving multidimensional boundary- value problems 

^ . by finite element method. An application of the algorithm to calculation of the 

Q ■ ground and first exited states of a Helium atom is given. 

O . Key words: Multistep adiabatic approximation, multidimensional boundary- 

value problems 

> 

Q\ ■ 1 Motivation 

oo ■ 

^^ ■ The adiabatic approximation is well-known method for effective study of few- 

^N I body systems in molecular, atomic and nuclear physics. On the base of pio- 

l/^ ■ neering work of Born and Oppenheimer [6] the method was applied in various 

^D I problems of physics, using the idea of separation of "fast" Xf and "slow" Xs 

variables [5 in Hamiltonian composed by fast and slow subsystems II{xf^Xs) = 
Hf{xf;Xs) + IIs{xs) with characterized frequencies Uf > cj^, for example in 
Henon-Heiles model [16 or quantum dot (QD) models [12]. 

Purpose of this paper is to present algorithm for generalization of the stan- 
;J] ■ dard adiabatic ansatz [15, 4 , 

{xf,Xs\nk) := Y^ , (f/K+i,f,)(f5,<+i|nfe), (1) 

for the case of multi-channel wave function when all variables treated dynam- 
ically p^ and to give a general scheme and examples of its application to cal- 
culation of ground and exited states of Helium atom [HI [3] . 
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For this reason we are introducing the step-by-step averaging methods in 
order to ehminate consequently ordered independent variables {x = {xf^Xs} = 
{xn y xn-i >^ ... )- xi}-^ G X = Xn U ... U Xi of subspace of coordinate 
space X C R^) and to improve accuracy of calculations of the parametric 
basis functions and corresponded matrix elements, and to have possibility for 
reducing computer resources in multi-dimension case by using in perspective 
the MPI technology. 

We present a symbolic-numerical algorithm for reduction of multistep adia- 
batic equations, corresponding to the MultiStep Generalization of Kantorovich 
Method [14J named below as (MSGKM), for solving multidimensional boundary- 
value problems with solutions subject to corresponding boundary conditions [7J 

H^m -2En,^m =0. (2) 

Here the Hamiltonian H = X]i=i H^+i-i of a quantum system presented by 
sum of parametric Hamiltonians Hi = Hi{xi;Xi-i^ ...^xi) of subsystems, de- 
pending on subset of independent variable Xi and parameters x^-i, ...,xi, and 
solutions satisfy to orthogonality and normalizing conditions 



{n[\ni) = / dXN'''dXilljl,{x)llJni{x)=Sn'ni' (3) 

For solving of problem (ISl)-© we propose multistep generalization of the stan- 
dard adiabatic ansatz ([1]) in the following form: 

Optimization of a convergence rate of the method is possible at an appropriate 
choice of characterized frequencies uon > ^n-i > ••• > oji of subsystems. 

2 Algorithm MSGKM 

Below we present symbolic algorithm MSGKM for generation of the boundary- 
value problems realizing multistep adiabatic expansion (|4]) in solving eigenvalue 
problem ([2]) in a symbolic form by using a Maple system. The examples of 
different versions of the algorithm are given in the next sections. 

Algorithm MSGKM 

Input: 

H = Xli=i ^A^+i-i is main Hamiltonian dependent on ordered variables x = 
{xn y xn-1 >-...>- xi}^ decomposed to sum of parametric Hamiltonians Hi = 
Hi{xi] Xi-i^ ..., xi), dependent on subset independent variable Xi and parameters 

Htljn^ -2En^^l)n^ =0, 



is main eigenvalue problem for calculation of unknowns 

i^m = \ni) ^ {x\ni) = i^mix) and 2En^ = £n^. 

Output: 

A set of Eq(/c), /c = 1, ..., A/", is a set of auxiliary parametric eigenvalue problems 
for calculation of V^^fc = V^it (^a^^ •••^ ^fc; ^^/c-i, ..., xi) and elfj = £lf;(x/e_i...xi), 
where ipm = V^ni and 2En^ = Sni are solutions of the main eigenvalue problem. 



Local: 

V^nt = V^it (^Ar,...,x/e;x/e_i,...,xi) and En^ = £nj = sl^J {xk-i...Xi) are solu- 
tions of the auxiliary parametric eigenvalue problems: 



Jk) 



Sk)( 



ik) 



Sk)j(k) 



(KM = f: 



0, 



dXN'-'dXN+l-ki^l^-^ i^nJ=Snlnk'i 



K+ll^fe) 

as: 



XnU...UXn 

X / {xk't Xk-ij ...jXi) are auxiliary parametric solutions defined 



the square brackets [, ] means a commutator: 

K+i| Hk.n'j^^^) =(n/e+i|i:ffe<+i)-K+iK+i)i:^fe. 



Eq(7V):={i7n,|n^) - £n,|n^) = 0, {rnV'Wn') = ^n^wj 



Eq(7V) ^{|n7v),£n^} 
for k:=N-l:l step -1 

Eq(fe):={(i7fc + £it+^) -sit^)K+i|nfc) 

+ En;^, K+l I [^fc, 4+l)J K+l |n/e) =0}. 

5: Eq(/c)^{«+,|n,),£it^} 
end for 



V^m = l^i), 2£:ni =^ 



- .(1) 



3 Statement of the problem for a Helium atom 

The Schrodinger equation for a Helium atom with total zero-angular momentum 
in hyperspherical coordinates [3]: = x^ ^ X3 = [0,7r],Qf = X2 G X2 = 

[0,7r],i? = xi G Xi = [0,+(X)), f = {X3 >- X2 >- xiY G X = X3 UX2 UXi 
reads as, 



{H2,{x2,]X2,Xi) ^ H2{x2]Xi) ^ Hi{xi) - 2Ei)^f i{x2„X2,Xi) = 0. 



Here the Hamiltonians Hi of subsystems consistent of differential operators by 
independent variables and multiplication operators of the Coulomb potential 
energy Vi of the three interacted particles with charges Za = —1, ^5 = —1, 
Zc = 2, including appropriate choice of weight factors: 

4 

^3(^3; ^2,^1) = ^ . 2 i^3(^3;^2,^l), 

xf sm X2 
TT/ N Id. d - . , 

^^3(^3; ^2,^1) = : ^ — smX3T^ h l/3(^3;^2,^l), 

smxs oxs oxs 

y xisin^X2 ZaZb 



2 \/T^^smx2Cosxs^ 
4 ^ 

i^2(^2;^i) = — i^2(^2;^i), 



d 2 d 
- — sm X2^- 
sm" X2 0x2 0x2 

^1 I Z^Zc ZjjZc 



H2{x2',xi) = - .2 ^r~^^^^^2^r~ + ^ + "^2(^:^2; ^1), 



^2(^2;^!) 



2 Vsinf cosf 

i:fi(xi) = Hi{xi), 

H (x ) = - — — x^— - A 

xl dxi ^ dxi x\ ' 



Solutions of discrete spectrum satisfy to orthonormalization conditions 

- / sinxsdxssin^ X2dx2xldxi'^i{xs J X2 J xi)'^j{xs J X2 J xi) = Sij 
8 Jx 



/x 
and subject to the boundary conditions 

^d'^i{xs,X2,Xi) 5 

lim x^ = 0, lim Xi^/i{xs,X2jXi) = 0, (5) 

. 2 a^i(x3,X2,Xi) . 5^i(x3,X2,Xi) 

lim sm X2 = 0, lim smxa = 0. 

2:2^0, TT dX2 X3^0,7r dXs 

4 Algorithm 1. Example of the conventional 
Kantorovich method. 

We consider one-parametric boundary- value problem with respect to fast Xf = 
{xs,X2} independent variables 



sm X2 ^ 



( ^^2 + H2{x2;x,) - ^4f (xi))^lf (X3,X2;X1) = 0, (6) 



/ sin xsdxs sin^ X2dx2'^il\x3j X2; xi)'^fj {^^ = Si^j^, 

^XgUXa 



(7) 



and conventional one by independent variables x = {xs >- X2 ^ xi}, 

l4.H3{x3;X2,Xi) 4.H2{x2;Xi) , £y . . o 1:^(1) VxrCl)/ \ n /q\ 

\ xf sm X2 XI I 

- I smx^dx^sin^ X2dx2x\dxi'^f^{x3,X2,xi)'^y^^ (9) 

with boundary conditions following from (|5]). 

In Step 1 we find the required solution of the problem (|6j) in the series 
expansion over the Legendre polynomials P^^ (cosxa) for each values of xi: 

•max 
H 

<'\x3,x2;xi) = ^P,^(cosx3)xli(^2;xi). (10) 

Substituting expansion ([TOj) into equation (J6]) and projecting with account of 
orthonormalization conditions of Legendre polynomials, we arrive to the one- 
parametric problem for unknown vector eigenfunctions, Xj-^i2{^2]Xi)i and cor- 
responded eigenvalue (potential curve), E\ (xi). 



-^-^ — - — sin X2^ hl + V'2(a;2,a;i) H ^^ 

sm X2 0x2 0x2 sm X2 

-\E'^^\x^)\ xfX{x2:x,) + — i- ^ T>S(x2;xi)xg,(x2;xi) = 0, (11) 

^ Jl=l 

V>f^^{x2]Xi) = / sin xsdxsPi^{cos xs)Vs{xs, X2; xi)Pj^{cos Xs) 

with boundary conditions following from ([5]). 

Substituting expansion (p!Q]) into (O, we have orthonormation conditions 

Y^ / sin^X2(ix2xi^U^2;xi)x^-^U^2;xi) = (^iui. (12) 

This one-parametric problem is solved with help of the adaptation of K ANTBP 
program [7J, named here as KANTBP 3.0. 

In Step 2 we find the solution of the problem ([8]) in the series expansion 
over solutions of problem (J6j) solved in the Step 1, 

•max 

^[l\xs.X2.x^) = ^ ^[l\xs.X2;x,)xtl{xr). (13) 

Substituting expansion ([13]) into equation (|8]) and projecting with account of 
orthonormalization conditions ([7]) of parametric basis functions '^■J{xs^X2'',xi) 
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Figure 1: Calculated potential curves in step 1 of equation in step 2. Circles on 
right panel note avoiding crossing points. 

from (p^Qj) calculated in Step 1, we arrive to the problem for unknown vector 



M) 



(1) 



functions, Xiji-^i.^i)-) ^^^ corresponding eigenenergy, £^|^ , 



-Y,{i2\[Hi,32)\xfl{x,)=Q, 

J2 = l 



{i2\ 



Hl,J2) 



1 d 



„5 4l;0;l/ 



il;0;l/ 



d 



"^^^'^ - #a^"^"^-- (-1) - ^S^(-i)a^ 



(14) 



with boundary conditions following from (J5j). Substituting expansion (p!3|) into 
([9]), we have orthonormation conditions 



2 -j /> 

J2 = l •'Xl 



(15) 



In (p!4|) we have definitions of elements of matrix of effective potentials {h 
0,1): 



A^'^^^''{xi)= / smx3dx3 sin x^dxi 



a'i*,^f(x3,X2;a;i)a'-i*''^(x3,X2;xi) 



axj^ 



9° 



dx\^lJ{xz,X2\Xi) 



dx{' 
'^fXx^,X2\xi). (16) 



A parametric derivatives of eigenfunction ^^- (x3,X2;xi) in (p!6|) are calculated 



with help of KANTBP 3.0 program together with corresponding integrals, where 
integration by variable x^, perform analytically by using of orthonormalization 
conditions of Legendre polynomials. As an example, some potential curves and 
matrix elements of effective potential matrix are shown in Figs. [T]and[2l One can 
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Figure 2: Calculated matrix elements in step 1 of equation in step 2. 



Table 1: Ground state Isls energy E^ and first exited 
of Helium atom (in a.u.) versus number i^^^ of basis 
^max q£ |.]^g Legengre polynomials 



state ls2s energy E^ 
functions and number 



,-max 


Isls: i5"^^ = 12 


Isls: i5"^^ = 21 


Isls: i5"^^ = 28 


ls2s: i5"^'' = 28 


1 


-2.895 539 01 


-2.895 551 19 


-2.895 552 76 


-2.139 935 68 


2 


-2.898 631 39 


-2.898 643 21 


-2.898 644 74 


-2.141 664 33 


6 


-2.903 643 86 


-2.903 655 95 


-2.903 657 51 


-2.145 700 22 


10 


-2.903 702 68 


-2.903 714 86 


-2.903 716 36 


-2.145 915 09 


15 


-2.903 708 49 


-2.903 720 68 


-2.903 722 17 


-2.145 957 35 


21 


-2.903 709 31 


-2.903 721 50 


-2.903 722 994 


-2.145 968 77 


28 


-2.903 709 31 




-2.903 722 997 


-2.145 970 28 


im 






-2.903 722 998 


-2.145 956 975 









-2.903 724 377 


-2.145 974 046 



see that these matrix elements have smooth behavior with respect to parameter 
x\ that achieve by imposing conditions of continuity of eigenfunctions with 
respect to parameter x\ in avoiding crossing points of potential curves E\^' {x\)^ 
shown by circles in Fig. (TJ where change number of zeros of corresponding pair 
of eigenfunctions by each of independent variables X2 and xz occurs after passing 
these points, as discussed in [3[|2]. As we show in Table [TJ such continuation via 
these points leads to increase of the convergence rate of expansion (p!3|) of total 
solution in calculation of ground and first exited states energy E^ with respect 
to previous one [3 . One can see from the Table [1] that convergence start from 
^max _ 2x is slow with rcspcct to upper variational estimation [9 . So, to improve 
convergence of calculation of the parametric basis functions from expansion (p!Q|) 
by number if^^ > 28, we begin to study in the next section the step-by-step 
averaging method for realizing calculation with a more high accuracy with help 
of Algorithm 2. Meanwhile, our upper estimation at z^^^^ = 28 for first exited 
state is lowing than result of [11 . On Fig. [3] we show the radial eigenfunctions of 
ground and first exited states. Note that, as following from asymptotic effective 
potentials (see Fig. [2]) the ground state solution has asymptotic in a vicinity 
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Figure 3: Radial eigenfunctions of ground and first exited states. 

triple-collision point R ^ including logarithmic terms that corresponding to 
Fock expansion [2j |13] . 



5 Algorithm 2. Example of MultiStep Gener- 
alization of Kantorovich Method (MSGKM) 

We examine a split sequence of boundary- value problems consists of the two- 
parametric problem by one of fast independent variables, Xf = {xa}, 



{H3{x3;x2,xi) - -E^^\x2,xi))^[l\x3;x2,xi) =0, 



(17) 



sinxsdxs'^l Uxs; X2, xi)"^) \xs; X2, xi) = S, 



Xa 



*3J3 5 



the one-parametric problem (|6])-([7]) by fast independent variables Xf = {xs >- 
X2} and conventional problem (|8])-([9]) by independent variables x = {^3 y X2 >- 
Xi} with corresponding boundary conditions following from (jSj. 

In Step 1 the two-parametric problem (pT|) with boundary conditions fol- 
lowing from (J5j) is solved for each values of xi and X2 with help of the adapta- 
tion of ODPEVP program [8 , named here as ODPEVP 2.0. The eigenvalues 



.(3)/ 



= ^(3) 



(potentials surfaces) El^\x2;xi) = El^\x2^xi) and their parametric deriva- 
tives are presented on Fig. ID One can see from Fig. HI the potential sur- 
faces are symmetric with respect to axis X2 = 7r/2, then the partial derivative 
dEl^\x2,xi)/dx2 = for is = 1,2, .... 

In Step 2 we find the solution of the one-parametric problem (J6]) in the 
series expansion over solutions of the problem (p!7|) solved in the Step 1: 



^lf{xs,X2;xi) = Y^ ^I\x3;x2,xi)x[li{x2;xi). 

i3 = l 



(18) 



Substituting expansion (p!8|) into equation (J6]) and projecting with account of or- 
thonormalization conditions (pT|) of parametric basis functions from Step 1, we 




Figure 4: a. The eigenvalues (potentials surfaces) E^ (^27^1)7 ^3 = 17 2,3,4 
of the problem ([T7j). b-c. their parametric derivatives dElJ {x2^xi)/dxi 
and dElJ {x2^xi) / dx2^ is = 1,2,3,4. d. their mixed parametric derivative 

dE[ \x2->xi)/dxidx2 



and eigenvalues (potentials curves) ElJ {xi): 



arrive to the one-parametric problem for unknown vector functions Xi.c.^ (^2; ^1) 



1 a . 2 5 



sin x-2 0x2 dx2 



eI'\x2;x^) 



2 sin X2 



(19) 



-i^£')(xi))xi^U^2;Xi)+^(z3|[i^2,j3)l^ 
/• I rr • \1 / /,2;10;10/ n 1 ^ • 2 42;00;10/ n 

(^al H2,j3) = I A3J3 (^2;^i) - — 2— ^::7sm X2A.'^.^' (X2;xi) 



sin^ X2 9^2 



^2;00;10/ n ^ \ 

with boundary conditions following from ([5j). Substituting expansion ([T8|) into 
([7j), we have orthonormalization conditions 



•iricix 

Y^ / Sin^X2(iX2X^3i2(^2;^l)Xj3J2(^2;^l) = ^i2j2- 



(20) 



J3 = l'^^2 



In (p!9|) we have definitions of elements of matrix of effective potentials: 



Alf^^--''-^'-^X2;x^)= f sinxsdx: 



^h+hxui^) 



5^2+^^^^^"^(x3;X2,Xi)S^2+^^^3^(x3;X2,Xi) 



(3). 



dx[^dx[' 



dx^2^^V 




Figure 5: Calculated elements A{2 ' (^2;^i) and A{^ ' {x2'-,xi) of matrix of 
effective potentials of Eqs. (p!9|) . 




Figure 6: Calculated elements A^^ ' (x2;xi) and A22 ' (2:2; ^1) of matrix of 
effective potentials of Eqs. (fT9|) . 

^J^*g)(x3;x2,xi) ^ *(3^)(x3;x2,xi). (21) 

Note that this problem is similar to problem (p!4|) from Algorithm 1, but ele- 
ments of matrix of effective potentials are calculated here with eigenfunctions 
(pTj) and their derivatives by parameters X2, xi by program ODPEVP 2.0 with 
accuracy 0{h'^^), where p is degree of approximation in a finite element grid [8 . 
As an example, some elements of effective potential matrix are shown in Figs. [5l 
[6] and El The eigenvalues (potential curves) El^\xi) of (fTO]) calculated here by 
program KANTBP 3.0 look as the same as in Figs. [T] calculated by Algorithm 
1. 

In Step 3 we find the solution of the problem (J8]) in the series expansion 
over solutions of problem (|6j)-([7j) solved in the Step 2: 



^\l\xs,X2,x,) = J2 ^\l\xs,X2;x,)x^X(^i). 



(22) 



J2 = l 



Substituting expansion (|22p into equation (jH]) and projecting with account of 
orthonormalization conditions ^ of parametric basis functions from Step 2, 



10 




3^"^ 3 

Figure 7: Calculated elements 74;l2 ' (^2; ^i) and A^^g ' (x2;xi) of matrix of 
effective potentials of Eqs. (p!9|) . 



(1) 



we arrive to the problem for unknown vector functions Xi^iii^i)- 



1 d 



xl dxi "*" dxi 



d . 2EfJ{x^ 



^(1) ,.(1) 



2E:;Mxi:i,(xi) (23) 



(i2| H,,n) ={a^^(x^ 



i2=i 
^1^^2.72 va^ij ^uj^ y^i) , 






into 



with boundary conditions following from ([5]). Substituting expansion 
(J9]), we have required orthonormation conditions 

•max 

If 

In (|23|) we have definitions of elements of matrix of effective potentials: 



(24) 



A^risci) 



sin X3 (ix3 sin X2 dx2 - 



X3UX2 



dx[^ dx\^ 

aQ^lf(x3,X2;xi) _ (2) 



dxi 



^•, (^3,^2;^i). 



Substituting expansion (p!8|) , we reduce matrix elements A^' -^'^^(xi) to integrals 



*2J2 



by variable X2 only via matrix elements A^^^^ ^'^^^^(x2; xi) calculated with help 
of improved parametric basis functions (p!8|) from Step 2: 



A, 



i;^i;ri 



h ri 



«3,J3 ki=0 kr=0 



h\ 



n\ 



h\{li-h)\kr\{n-kr)\ 



(25) 



• 2 J ^2;0fei;0fe,./ \ 

sm X2dx2A^'," '^{X2;xi) 



d''-'''xtli^2;xi)d-^-'^^xfXix2;xi) 



dx['-''' 



dx[ 



ri—kr 
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Note that set of Eqs. (|23|) is similar to Eqs. ([H|) from Algorithm 1, and the 
potential curves ElJ (xi) and elements of matrices A-^j^^^{xi) calculated by 
Algorithm 2 look as the same as in Figs. [1] and [2] calculated by Algorithm 1. 
We can wait that using expansion (p!8|) over the two-parametric basis functions 
(pT|) will have a higher rate of convergence and give corresponding elements of 
matrices of effective potentials (|25]) with better accuracy in comparison with 
expansion (p^Qj) over the Legendre polynomials. However, matrix elements have 
a more comphcate structure and additional numerical integration by variable 
xs performed in Step 2 is needed in comparison with ([16]). 



6 Conclusions 

In this paper we presented a symbolic algorithm for reduction of multistep adia- 
batic equations, corresponding to the MultiStep Generalization of Kantorovich 
Method, for solving multidimensional boundary-value problems and consider 
examples of its application to Helium atom calculation. Achievement of this 
approach consist in facts that, on each step solution subject to boundary condi- 
tions, elements of matrix of effective potentials calculated with controllable ac- 
curacy have smooth behavior with respect to parameters like that in spheroidal 
coordinates [1 . These facts together with consistency of Kantorovich expansion 
in a vicinity triple-collision point R ^ including logarithmic terms corre- 
sponding to Fock expansion [2l [13] provide as shown above a reasonable rate of 
convergence of these expansions and upper estimations of energy eigenvalues [16] 
Moreover, asymptotics of these expansions at large value i? ^ oo in vicinities of 
pair collision points of limit of separated atom are compatible with asymptotic 
states needed for solving a scattering problem below three body threshold, as 
shown in papers [2", 13^ . 

Elaboration final version of Program KANTBP 3.0 for solving the problem 
with respect to unknowns (i.e. calculation of improved parametric basis func- 
tions in Algorithm 2) from Steps 2— (n-1), with matrices of variable coefficients 
calculated and presented above is in progress. 

Generalization of MultiStep Kantorovich method presented above reduce to 
the set of 2N — 1 of multiparametric eigenvalue problem for set of ~ 10 ordinary 
second-order differential equations that can solve naturally by each of A^ — 1, 
A" — 2, ..., 1, independent parameter using MPI and/or GRID technology that 
will be elaborated in our further investigations. 

The computational scheme, the SNA, and the complex of programs allow 
extension for the analysis of spectral characteristics of both electron(hole), im- 
purity and excitonic states in nanoscale quantum-dimensional models like QWs 
[T7], QWrs[T8]. QDs [12] with different geometry of structure and spatial form 
of confining potential and external fields. 
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